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Abstract 

This Letter presents an investigation on the effects of mutual coupling 
in a metamaterial comprising two sets of electric-LC (ELC) resonators 
with different resonance frequencies. Through simulation and experiment, 
it is found that the two resonances experience significant shifting and 
weakening as they become spectrally close. An equivalent circuit model 
suggests that inductive coupling among the two resonator sets is a primary 
cause of the change in the resonance properties. This study is fundamen- 
tal to designing metamaterials with an extended bandwidth or spatially 
variable response. 



1 Introduction 

A metamaterial defines a group of resonators that collectively exhibit a strong 
electric and/or magnetic resonance. These resonators operate in the effective 
medium regime owing to their subwavelength dimensions. Thus, their response 
can be characterized by a homogenized permittivity and permeability, which 
can be controlled via the shape, size, and materials of the resonators. In the 
last decade, metamaterials have appeared in various applications across the 
electromagnetic spectrum [T]. 

To achieve desirable properties, the knowledge on inter-resonator coupling is 
significant, since this near- field behavior directly influences the effective response 
of metamaterials. Earlier studies have focused on coupling among identical res- 
onators [2J [3J 21 O [5] ■ Recently proposed metamaterials contain different unit 
cells to overcome the inherent narrow-band limitation [71 [51 [HI [TU1 EJ [T2J [T3] . 
Additionally, transformation optics requires a spatial change in the metamate- 
rial response [Ml [15] . Hence, the study on the interaction of different resonators 
become vital. 
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In this work, mutual coupling in a dual-resonance metamaterial, composed 
of different electric-LC (ELC) resonators [16j E3 HB] > is explored. The coupling- 
effect is observed through simulation and experiment. A microwave metamate- 
rial is used for ease of continuous resonance tuning via integrated varactors [19] . 
A double-tuned resonance circuit along with a transmission line model provides 
further insight into the coupling phenomenon. 



2 Experiment 

A unit cell of the varactor-integrated ELC resonator is illustrated in Fig. UJa). 
The dimensions of the resonator are as follows: a — 13 mm, d = 11 mm, 
I = 1.8 mm, w — 1.7 mm, b — 0.8 mm, g = 0.4 mm, and s = 0.3 mm. Each 
resonator is made of 35-^m-thick bare copper, supported by an epoxy FR4 
substrate with a thickness of 0.8 mm, a dielectric constant of 4.5, and a loss 
tangent of 0.02. The entire metamaterial panel comprises 20x20 resonators. 
The varactor is Infineon BB837 with a tunable range from 9.5 pF down to 
0.58 pF for the reverse bias from to 24 V. The varactors in each row are 
connected in parallel to an independent dc source. Fig.QJb) shows the fabricated 
metamaterial including dc feeds. 

The experiment is carried out in an anechoic chamber with a vector network 
analyzer connecting to transmitting and receiving horn antennas. The sam- 
ple transmission is measured and normalized by the free-space transmission. 
The simulation is performed using Ansoft HFSS employing the finite-element 
method. Periodic boundary conditions are used for the transverse boundaries 
to replicate an infinite 2D resonator array. The varactor is substituted by a 
lumped capacitor, whose capacitance is available from the SPICE model. 

The experiment comprises three cases. First, all of the varactors in the array 
are connected to a dc bias that changes from to -24 V with an interval of - 
4 V (case S), the results of which are elaborated in [J5]. Next, the varactors 




Figure 1: (a) ELC resonator with a varactor loaded at the central gap. The dc 
bias lines are positioned at the four corners. The dotted lines bound a unit cell, 
(b) Fabricated planar metamaterial containing 20x20 resonators. The biasing 
strips are visible on the left and right. 
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Figure 2: Transmission magnitudes from experiment (a,b), simulation (c,d), and 
circuit model (e,f). (a,c,e) The odd-row bias is V, whilst the even- row bias 
varies from to -24 V (case DO). (b,d,f) The odd-row bias is -24 V, whilst the 
even-row bias varies from to -24 V (case Dl). 




Figure 3: Simulated resonance frequencies as a function of bias voltage for three 
different bias cases: S, DO, and Dl. 



in the odd rows are grounded, whilst those in the even rows are connected to 
a bias varying from to -24 V with a step of -4 V (case DO). Finally, the 
odd-row varactors are pulled to -24 V, and the even-row bias varies from to 
-24 V (case Dl). The experimental results for cases DO and Dl are shown in 
Fig. (2ja,b), whilst Fig. [2](c,d) are the corresponding simulation results. The 
simulation and experiment yield similar results with small discrepancies. The 
transmission magnitudes for both cases DO and Dl reveal the lower and higher 
resonances, caused by lower and higher absolute biases, respectively. 

The resonance frequencies as a function of the bias voltage are depicted in 
Fig. [31 In both cases DO and Dl, the lower and higher resonances experience 
a slight blueshift, relative to the resonance of the single-biased metamaterial 
(case S) at the same bias voltage. This blueshift is related to the fact that 
the inter-cell coupling, typically reinforced among the ELC resonators with the 
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same resonance frequency, is weakened as half of the resonators are spectrally 
tuned away. 

In both cases DO and Dl, the strength and linewidth of the higher resonance 
are noticeably smaller than those of the lower counterpart, and the difference 
is more pronounced when the two resonances come closer to each other. This 
effect, also observable in |S], is attributed to the different oscillation phases 
in the lower-biased resonators. At the lower resonance frequency, the in-phase 
oscillation induces constructive current in the higher-biased resonators. Slightly 
above the lower resonance frequency, the out-of-phase oscillation [5Q1 123] induces 
destructive current in the higher-biased resonators. Hence, the imbalance in the 
two resonances is observed. 



3 Equivalent circuit model 

A double-tuned resonant circuit model shown in Figgis utilized to gain further 
insight on the coupling mechanism of the complex dual-resonance metamaterial 
[22] . It is composed of two inductively coupled series RLC circuits. This sim- 
plification lies on the basis that an ELC resonator can be approximated by a 
series RLC circuit with a Lorentzian response and that the magnetic flux link- 
age exists between neighboring resonators. The two resonant circuits in Fig. |4] 
share the same ac source, since their equivalent ELC resonators are uniformly 
excited. 

Through Kirchhoff 's current law, the system of equations for the circuit in 
Fig. S] is given as 
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with complex impedances Z{ 12 } = ^{1,2} + J w -k{i.2} + and Zm = 

juiM . The sign for the mutual inductance, M , depends on the magnetic field 
activity. It is positive if the overall mutual induction is constructive, and be- 
comes negative otherwise. From Eq. Q] the impedance seen by the source Vq 
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Figure 4: Double-tuned resonant circuit. The phasor voltage and current are 
represented by Vq and I{i,2}, respectively. 



4 



can be described as 



Z c = 



V 



ZiZ 2 — 7i\ 



(2) 



Ii +I 2 



Zi + Z 2 — 2Zm 



Transmission line theory is used in conjunction with this circuit model to de- 
termine the transmission characteristics [221 123] ■ The substrate introduces a 
parallel impedance Z s to Z c , and also causes a change in the transmission mag- 
nitude at its back surface. Hence, the model complex transmission coefficient is 
given as 



where the free-space and substrate intrinsic impedances, Zq and Z s , equal 377 
and 177.72 fi, respectively. 

Eq. [3] contains seven variables: R{l.2}^ C{i.2}i £{1,2}? and M. Here, R\ = 
i?2 = Rs and Li = L2 = L s , because of an identical resonator structure. The 
capacitances, Cn 2}i equal the variable capacitances C{di,d2} m series with the 
structure capacitance C' s . Although C{di,d2} ar e correlated with the varactor 
capacitances, and C s , and L s can be readily obtained [19] . using these values do 
not yield satisfying results. This is because the simplified circuit model repre- 
sents a large resonator array, and thus the individual capacitors and inductors 
in the model are different from those in the array. Hence, only a relative change 
is deduced from these known values. 

Based on these assumptions, the fitted values are as follows: R s — 12 fi, 
C s = 92 fF, L s = 36.5 nH, M = 2.8 nH. The variable capacitance C d = {1.9, 
0.62, 0.32, 0.21, 0.16, 0.13, 0.12} pF for the reverse bias varying from to 
24 V in 4-V steps. It is clear from the transmission results in Fig. [2je,f ) that 
the circuit model can mimic the behavior of the dual-resonance metamaterial. 
The same parameter set also produces an agreeable transmission magnitude for 
the single-resonance metamaterial. A small difference in the resonance shape is 
mainly caused by the high-frequency dipole resonance |23j . which is not taken 
into account by the model. 

The mutual inductance M implies the magnetic flux linkage between the two 
resonator groups. By changing this value, the influence of magnetic coupling on 
the metamaterial response can be observed. Fig. [5] depicts the model transmis- 
sion magnitudes following a change in M . As M is reduced, the strength of the 
lower and higher resonances decreases and increases, respectively. Additionally, 
the lower resonance undergoes a blueshift, whilst the higher resonance under- 
goes a redshift. This effect can be ascribed to a reduction in the constructive 
and destructive magnetic coupling among the different resonators at the lower 
and higher resonance frequencies, respectively. At M = 0, where the two sys- 
tems are decoupled, the two resonances become comparable in strength. Hence, 
a difference in the resonance strengths is responsible purely by the inter-cell 
inductive coupling. 



T = 



2Z c ff 2Zq 



Z c ff = Z C ||Z. 



(3) 



Z c ff + z z s + z ' 
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Figure 5: Transmission magnitudes for various mutual inductance M obtained 
from the circuit model. The odd- and even-row biases are fixed to and 4 V, 
respectively 



4 Conclusion 

In conclusion, the effects of mutual coupling in a dual-resonance metamaterial 
have been studied through the simulation, experiment, and interpreted through 
a circuit model. As a result of constructive and destructive interferences in the 
magnetic activity, the coupled resonances are unbalanced and shifted consid- 
erably. This study is fundamental towards designing extended-bandwidth or 
spatially dependent metamaterials, particularly relevant to the emerging field 
of transformation optics. The results from this study are applicable to other 
types of resonators possessing a similar resonance mode in any frequency range. 
Future work involves minimising the coupling effect via optimal designs. 
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